function [act logact] = gmmactivBayesVerbeek(mix, data)


K = mix.ncentres;
N = size(data.Y, 1);

m = mix.varposterior.m;
b = mix.varposterior.b;
W = mix.varposterior.W;
v = mix.varposterior.v;
ElndetLambda = mix.varposterior.ElndetLambda;

act = zeros(N, K);
logact = zeros(N, K);
for k = 1:K
    [act(:, k) logact(:, k)] = gaussBayesVerbeek(m(k, :), b(k), W(:, :, k), v(k), ElndetLambda(k), ...
        data.Y, data.C);
end
